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ABSTRACT 



We study the dynamic spin correlation function of a spin one antiferromagnetic 
chain with easy-plane single-ion anisotropy. We use exact diagonalization by the 
Lanczos method for chains of lengths up to N=16 spins. We show that a single-mode 
approximation is an excellent description of the dynamical properties. A variational 
calculation allows us to clarify the nature of the excitations. The existence of a two- 
particle continuum near zero wavevector is clearly seen both in finite-size effects and 
in the dynamical structure factor. The recent neutron scattering experiments on 
the quasi-one-dimensional antiferromagnet NENP are fully explained by our results. 
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I. INTRODUCTION 



It was first argued by Haldane [1,2] that generic spin- S one-dimensional Heisen- 
berg antiferromagnets have an excitation gap for integer S. This picture is quite dif- 
ferent from the usual, higher-dimensional, picture of antiferromagnets with massless 
Goldstone magnons. On the theoretical side, there is now convincing evidence from 
numerical studies [3-10] of the S=l Heisenberg chain that it has a nonzero gap in 
the thermodynamic limit. Exact diagonalizations by the Lanczos method have been 
able to reach 18 spins and Monte-Carlo simulations extend the range to 32 spins. 
All these data are suggestive of a Haldane gap. This leads to a ground state spin 
correlation length which is finite and about six lattice spacings. Such a picture is in 
striking contrast with the 5= 1/2 solvable chain which is gap less and whose ground 
state has algebraically decaying correlations. A spin 1 Heisenberg model with bi- 
quadratic exchange has also been discovered with an exactly solvable ground state 
and a nonzero gap [11], thus reinforcing the belief in the Haldane conjecture. 

On the experimental side there are several candidates to exhibit this quantum 
gap. The first experimental evidence came from neutron-scattering on CsNiCls [12- 
15]. In this compound there are chains of spin one Ni 2+ ions with superexchange 
through Cl~ anions. There is however moderately small interchain couplings that 
complicate the picture: this is revealed in recent experiments on the related com- 
pound RbNiCl3 [16]. The best candidate so far seems to be Ni(C2H 8 N 2 )2N02C104 
(NENP) [17-20], as shown by inelastic neutron-scattering (INS) and magnetization 
measurements. The ratio of the interchain to intrachain magnetic coupling is es- 
timated to be J' I J ~ 4 x 10 -4 . No transition to Neel order is found down to 
1.2 K which is consistent with the hypothesis of a finite zero-temperature correla- 
tion length. Nickel ions have spin 1 and are described by the following anisotropic 
Heisenberg Hamiltonian [21]: 



Best fit of INS gap values leads to J/ks = 43. 5K and easy-plane D/J = 0.18 [10]. 
There is also evidence for a smaller in-plane anisotropy that can be described by 
adding to H a term E X]J(>S?) 2 ~~ O^T) ]• This perturbation is small and will be 
mostly ignored in the remainder of this paper. Its qualitative role will be discussed 
in section V. 

The first INS measurements [17-19] have concentrated on the neighborhood 
of Q=7r (where Q is the wavevector along the chain) and showed the existence of 
two gaps: one for the in-plane (IP) magnetic fluctuations and a higher one for 
out-of-plane (OP) fluctuations. This splitting is due to the presence of a sizable 
easy-plane single-ion anisotropy D. The dispersion of the magnetic excitations was 
then studied in the range Q/ir = 0.9 — 1.0. Recent experiments have extended our 
knowledge throughout the whole Brillouin zone [22]. Detailed theoretical work is 
required to test the hypothesis that the simple model Hamiltonian (1.1) is able to 
reproduce the dynamical properties of NENP. Two approaches have been followed 
up to now: small cluster numerical studies [3-10,23] and the use of effective field 
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theories [24,25]. The previous numerical studies were limited to spectral calculations 
or have ignored the presence of in-plane anisotropy. 

In this paper we present the results of a study of the dynamical structure factor 
S(Q, u>) of chains of length up to 16 spins by means of an exact diagonalization 
method. We also apply a variational technique proposed in ref.[26] to the anisotropic 
chain. This physically motivated method reproduces extremely well the numerical 
results. The results presented are valid in the absence of an external magnetic field 
and in the zero temperature limit, i.e. for temperatures well below the gap. In sect. 
II we present some general properties of the spin-1 chain with easy-plane single-ion 
anisotropy. In sect. Ill we explain the methods used to obtain dynamical quantities. 
Sect. IV contains our results. They are extremely well reproduced in terms of a 
physically appealing single-mode approximation. The results of our variational 
calculation are presented there and we use them to understand the nature of the 
elementary excitations near Q = and Q = it. We show in sect. V that they are in 
very good agreement with existing neutron data and we suggest an additional test 
of the theory. Sect. VI contains our conclusions. 

II. GENERAL PROPERTIES 

We briefly review some known results about the isotropic Heisenberg spin 1 
antiferromagnetic (AF) chain: 

N 

H = jJ2Si-S i+1 . (2.1) 

i=i 

The vectors Si are quantum spin operators satisfying the SU(2) rotation algebra 
with length Sf = 2. They are located along a one dimensional lattice of N sites 
with periodic boundary conditions. The exchange integral J is positive in the AF 
case. For any finite lattice the ground state is a singlet and the higher-lying levels 
have energy increasing with increasing spin as is known rigorously [27]. Above the 
singlet ground state one finds a triplet state presumably for all values of N. In 
an AF quantum magnet, in the case of broken symmetry one expects that in the 
thermodynamic limit the triplet becomes degenerate with the ground state as do 
other states with spin S=2,3,. . . , in order to form the degenerate ground state of 
the infinite volume limit. This can occur only for dimension greater than 2 (or 
equal). Haldane has argued that the spin-1 chain is quantum disordered and that 
the singlet-triplet gap remains nonzero in the thermodynamic limit. This should be 
true for any integer value of the spin. This argument is based on a mapping of the 
infinite-spin, semiclassical limit of the quantum chain onto a 0(3) nonlinear a model 
[1]. In the spin-1 case of interest, numerical studies clearly points towards a gap 
close to ~ 0.41J [3-10,23]. The physics of the 0(3) nonlinear a model is that of a 
triplet of massive bosons with nontrivial scattering properties. This suggests a very 
simple effective field theory: a free theory of three massive bosons [25]. Another 
route starting from an integrable model leads to an effective field theory of three 
massive fermions [24]. 
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In such approaches, one has to adjust the gap values which are no longer 
deducible from the microscopic model but one can obtain simple and definite pre- 
diction on dynamical quantities. However both are approximate theories and their 
respective domain of validity is difficult to assess. On the contrary, finite chains 
calculations offer unbiased theoretical predictions provided one is able to carefully 
control finite-size effects. Concerning the gap values this can be achieved by the use 
of the so-called Shanks transformation, suited to the removal of exponential tran- 
sients in sequences of finite-chain data [28] . It has been found that the convergence 
towards infinite-volume is very good [6,10,23] as expected since we are dealing with 
a massive theory. 

As shown with the Perron- Frobenius theorem [29], the lowest-lying triplet has 
wavevector Q = tt, while the singlet ground state has Q = 0. Most theoretical and 
experimental studies have concentrated on this Q = tt part of the spectrum. In 
the nonlinear a model, there are no bound states and, if we believe that it is the 
effective theory of the spin 1 chain, this implies that the Q — gap is due to two 
massive Q = tt particles and thus twice the Q = tt gap [2]. This has been checked 
by a quantum Monte-Carlo simulation in the isotropic case [8] and Lanczos studies 
[23] have shown that this property persists in the presence of anisotropy (as long 
as it is not too large) . As a consequence, the states near Q = should be members 
of a two-particle continuum contrary to the states near Q — tt that should appear 
as long-lived well-defined excitations. INS measurements for Q = [30] reveal a 
vanishing structure factor in this regime. The curve of the lowest excited triplet at 
wavevector Q is thus bell-shaped but asymmetrical with respect to Q = tt/2. 

Let us now discuss the influence of easy-plane single-ion anisotropy, i.e. adding 
to Hq a term of the following form: 



The full rotational symmetry of Hamiltonian (2.1) is broken to rotational symmetry 
around the z axis. Only the z component S z = S z of the total spin is conserved. 
There is also a discrete symmetry S z — > — S z that is preserved in the Hamiltonian 
(2.2). As a consequence of the anisotropy, the first excited triplet state is split into 
a higher-energy singlet S z = and a lower-lying doublet S z = ±1. These three 
states retain their wavevectors unchanged (Q = tt) with respect to the D = case 
(by continuity) since the Hamiltonian (2.2) still possesses translational symmetry. 
The Haldane gap is split into two components: one gap G^ between the ground 
state (with S z = 0) and the doublet S z = ±1 and one gap inside the S z = 
subspace between the first two levels. decreases while G^ increases with 

increasing D as shown quantitatively in ref.[10]. 

We expect this simple picture to be true for all wavevectors < Q < tt since 
the triplet states with arbitrary Q will be split also by the anisotropy. There is 
thus two distinct magnetic modes throughout the Brillouin zone: one IP (in-plane) 
mode (doubly degenerate) and one OP (out-of-plane) mode. The IP mode is seen 
in the S z = ±1 sector while the OP mode appears in the S z = sector. The 
dispersion of the IP mode has been obtained for various anisotropies including that 




(2.2) 
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of NENP in [23]. INS experiments [17-19,22] have resolved these two modes in the 
neighborhood of Q = tc: the OP mode is found at 2.5 meV and the IP mode is even 
further resolved into two components at 1.05 meV and 1.25 meV. This splitting is 
due to a small in-plane anisotropy of the type Ej^^Sf) 2 — (Sf) 2 ] that lifts the 
degeneracy of the doublet states S z = ±1. These values of the gaps (ignoring the 
in-plane anisotropy) leads to J = 44K and D/J = 0.18 [10]. Detailed interpretation 
of INS requires a calculation of the structure factor S(Q,uj). 

III. EVALUATING DYNAMICAL QUANTITIES 

We describe in this section the methods used to compute the dynamical struc- 
ture factor: 

S aa (Q,uj) = Jdt e*"* < 0|S? g (t)£§(0)|0 >, a = x,y,z. (3.1) 

Here we denote the ground state by |0 > and the components of the Fourier trans- 
form of the spin vectors are £g(t) in the Heisenberg representation. We first write 
the structure factor as: 

S aa (Q,u) = Y / S("-(en-e )) |<n|^|0>| 2 , (3.2) 

n 

where the sum over \n > means over all excited eigenstates of the system and e n 
denotes the energy of \n >. In a finite system calculation, one obtains a finite set 
of delta functions with weights given by the matrix elements appearing in Eq.(3.2). 
The Lanczos method which is commonly used to obtain the ground state wavefunc- 
tion and the first few excited levels is suited to the study of the dynamical properties 
[31]. 

We proceed as follows: 

i) One needs to know first the vector |0 >. We use the Lanczos algorithm applied 
to the Hamiltonian H. In fact, any kind of algorithm can be used at this level. 

ii) One constructs the state: 

|$ >=Sq|0>. (3.3) 

This state |$o > is then used to build a new Lanczos sequence of states |<fr n > by 
the standard procedure of applying the Hamiltonian to \<&k > and orthonormal- 
ization with respect to the last two vectors \<&k > and |$fc-i >. The Hamiltonian 
H is then tridiagonal in the basis {|$ n >}• The nonzero elements of H in this 
basis are directly provided by the orthonormalization procedure. We exhaust the 
corresponding subspace and store the tridiagonal form of the Hamiltonian. 

iii) One diagonalizes the tridiagonal Hamiltonian by a standard routine. This leads 
to a set of eigenenergies e n as well as the eigenvectors whose first coordinate in the 
Lanczos basis are precisely the overlap matrix elements < n|$o > that appear in the 
definition of the structure factor Eq.(3.2). We thus obtain the weight of each delta 
function peak. This has to be done separately for each Q as well as for a = x, z. 

In the isotropic D=0 case the ground state is a singlet: Sq =0 |0 >= and thus 
S aa (Q = 0, to) = 0. The vanishing with Q will occur quadratically. In the pres- 
ence of easy-plane anisotropy the ground state is invariant only under z rotations: 
S£ =o |0 >= and thus S ZZ (Q = 0,u) = 0, while S xx (= S yy ) will be nonzero at 

Q = and of order 0(D 2 ). With in-plane anisotropy EJ2 t [(S?) 2 - (Sf) 2 ] even 
S ZZ (Q = 0) will be nonzero and 0(E 2 ). 
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IV. RESULTS FOR S(Q,u) 

We have computed the structure factors S xx (Q,ui) and S zz (Q,uj) for chains 
of lengths N=4,6,8,10,12,14,16. Some Q = and Q = it parts of the spectrum 
have been computed also for N=18. We have concentrated on the case D/J — 0.18 
quantitatively relevant to the study of NENP. Our results extend smoothly for 
not too large anisotropy D < J. For D = we reproduce previous findings [8] 
concerning the lowest excited levels. Previous Monte-Carlo measurements of the 
dynamical structure factor [9] for D = are also compatible with our data. 

According to the definition (3.2), we will discuss first the peak positions and 
then the corresponding weights. What we observe is that for all fixed wavevector 
Q, as a function of u> the lowest lying peak concentrates almost all of the spectral 
weight. This remark will be made quantitatively precise at the end of this section. 
This lowest-lying peak occurs when the frequency 10 matches the energy of the first 
excited level with the right quantum numbers. 

In the correlation S xx we find the excitation spectrum from S z = to S z = 
±1 that we obtained previously [23]. This is plotted in fig.l as a function of the 
momenta by crosses of various sizes. For Q = and Q = tt there is enough 
data to allow an extrapolation to the thermodynamic limit by use of the Shanks 
tranformation. We have also computed the these gaps (Q = and Q = tt) for a 
N=18 chain. The extrapolation leads to 0.301J at Q = tt and 0.986J at Q = 0. In 
the interval [0, 7r] there is only the middle of the Brillouin zone Q = n/2 where we 
obtain several data points from N=4,8,12,16. Shanks extrapolation leads then to 
2.75J for Q = n/2. 

In the correlation S zz we find the excitation spectrum inside the S z = sub- 
space. The lowest lying excitation is plotted in fig.l with empty octogons. At Q = n 
the finite-lattice data extrapolate to 0.655 J. The convergence in this case Q = n is 
very good as expected from a massive theory yielding well-separated eigenvalues in 
the finite systems [23] . On the contrary at Q = the convergence is extremely bad 
and one can only suggest a value of 0.60 ± 0.05 J for this gap. This is clearly due 
to the presence of a continuum of states starting immediately above this gap. Such 
a phenomenon is expected if one describes the spectrum near Q = as due to two 
particles with momenta near Q = n [8]. 

We note that the excitations in this sector (zz) are above the in-plane excita- 
tions in the neighborhood of Q = tt but they cross for Q/tt ~ 0.75 and on the top of 
the dispersion curve it is the in-plane mode which has highest energy: 2.75 J against 
2.65J (extrapolated value) for the out-of-plane mode. This behaviour persists in 
the region 0<<5<7r/2of the Brillouin zone. 

To gain understanding of the nature of the excitations we have performed a 
variational calculation along the lines proposed by Gomez-Santos [26] . One discards 
the states with parallel spins, either if they are nearest-neighbors or separated by 
any number of zeros. The typical states in this subspace are of the form: 

| . . . tit 1 00 ti - • - > (4-1) 

There is strict Neel ordering of the S z = ±1 sites but there can be any number of 
intermediate zeros. In this subspace the true degrees of freedom are then the " spin- 
zero defects" i.e. sites with Sf = 0. These domain walls are then represented as 
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fermions (which is a natural way of enforcing no double occupancy). The fermionic 
Hamiltonian can then be treated by approximate methods. It has been shown that 
a simple Hartree-Fock decoupling leads to a good approximation of the spectrum 
for D = 0. The best results are obtained by using the variational improvement 
introduced by Gomez-Santos where one allows a small admixture of states with 
parallel spins. The Haldane gap, as computed by this method, is found to be 0.45J 
for D=0. It is straightforward to include anisotropy of the form D^^Sf) 2 . We 
have obtained a spectrum of massive fermions E(Q) that is plotted as a continuous 
curve in fig.l. This curve should be the excitation spectrum in the S z = subspace 
due to the restriction to the subspace (4.1). The gap at Q = it is found to be ~ 0.70 J 
when D/J=0.18. The single-fermion curve E(Q) reproduces extremely well the OP 
dispersion in the region n/2 < Q < n. The IP dispersion will be discussed later on. 



Figure 1 

Dispersion of magnetic excitations S z = and S z = ±1 at D = 0.18 J for momenta 
ranging from to it. The continuous line is the result of a variational calculation. 
The dashed line is the corresponding edge of two-particle excitations. 
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We now turn to the discussion of the spectral weight associated to the peak 
positions. For all chain lengths studied for all wavevectors in the interval [~ 0.3 — 
1.0] x 7r we find that both S xx and S zz are dominated by a single peak as a function 
of the frequency. This peak concentrates at least 98 per cent of the spectral weight. 
Multimagnon contributions are thus extremely small in the whole range [ps 0.3 — 
1.0] x 7r. As typical examples we plot in fig. 2 S XX (Q = tv,uj)+S zz (Q = 7r, uj) 
from the N=16 chain. The peaks have been broadened for clarity. In fig. 3 we 
plot the Q = 7r/2 case showing the interchange of the two modes. A single-mode 
approximation will be an extremely good description of the physics in this part of 
the Brillouin zone. The simple picture described above breaks down for small values 
of Q. In the N=16 chain, for Q=37r/8 there is still the one-peak structure while for 
Q=7r/4 there are several peaks as a function of uj: see fig. 4. In the N=14 chain we 
find that the continuum appears already for Q=27r/7 but is not there for 37r/7. It is 
difficult to give a precise value of the wavevector at which this phenomenon appears 
since the discretization of the momenta imposed by the chain length is coarse. As 
discussed below, this can be interpreted as evidence for a two-particle continuum. 
Our present data suggest that the continuum sets in at ~ 0.37T for both IP and OP 
sectors. We see no reason why the continuum boundary should the same for both 
modes although we are unable to see any quantitative difference in our present data 
between IP and OP modes in this respect. 

Let us discuss the nature of the two-particle excitations. If a single-particle ex- 
citation has a dispersion E(Q), the lower edge of the continuum of two-particle exci- 
tations is given by E^ 2 \Q) = m\rvK{E{K) + E(Q — K)) in the non-interacting case. 
Such a free picture seems to apply to the spin 1 chain [1,2,24]. In the anisotropic 
system it is possible to construct several continua since there are two distinct modes 
(IP and OP). 

In the S z = ±1 subspace two-particle states are obtained from one S z = ±1 
excitation and one S z = excitation. The gap of this continuum at Q = will be 
the sum of the IP and the OP gap at Q = it i.e. ps 0.301J+0.655 J in agreement with 
our extrapolation of 0.986J. The continuum is also seen in the poor convergence of 
finite-size data. Its progressive build-up is seen in fig. 4. 

In the S z = subspace there are two ways of building two-particle states: 
either with two S z = states or with one excitation S z = +1 and one excitation 
S z = — 1. In the first case we can obtain an approximation for the continuum 
boundary by using the fermionic method cited above: the curve E (2 \Q) for two 
fermionic excitations belonging to the S z = sector is plotted as a dash-dotted line 
in fig.l. The two-fermion continuum is clearly above the OP (and IP) mode for 
< Q < 7r/2. However the second possibility which is the continuum (S z = +1) 
+ (S z = —1) has a Q=0 gap that is twice the gap at Q=7r for the IP mode i.e. 
w 0.6 J. This agrees with our estimate ps 0.6 J for the gap in the OP sector at Q=0. 
This continuum is thus the lowest-lying one. This reasoning shows also that the 
OP (S z = 0) mode should be below the IP mode for Q ps since the IP gap is then 
given by ~ 0.301J + 0.655 J. In fact, as stressed above, we have good evidence that 
the crossing of the two modes occurs quite near Q/n ps 0.75. 
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Figure 2,3,4 

The structure factor S^ XX \Q, uj) + S^ ZZ \Q, uj), for Q = n, n/2, n/4 respectively, as 
a function of the frequency uj (in units of J) for chain length N=16. In fig. 2 each 
of the two functions 5'( aQ ) has a single-peak structure: a single-mode approximation 
works very well. The IP mode is doubly degenerate and appears only in while 
the OP appears in S^ zz \ln fig. 3 the IP and OP modes are interchanged with respect 
to fig. 2. In fig. 4 a single-mode approximation no longer holds and the two-particle 
continuum is seen: spectral weight is transferred to multimagnon excitations. 
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The agreement between the ab initio data and the variational calculation con- 
firms the relevance of domain walls in the whole Haldane phase, as is also pointed 
out in [32]. The disordered-flat phases obtained in the solid-on-solid picture of the 
spin-1 chain [33] also correspond to the subspace of spin-zero defects which is so 
successfull in the description of the elementary excitations. 

Finally we quote results for the static structure factors S^ aa \Q) obtained by 
integration over frequency. In fig. 5 these quantities are plotted in a logarithmic 
scale as a function of momentum. Due to our limited chain length, we cannot make 
any firm statement about the prediction [1] of square-root Lorentzian behaviour 
near Q = n. Our data are compatible with previous studies [34,35] concluding to 
a correlation length « 6 lattice spacings in the D = case: the correlation length 
does not change much when D/J = 0.18. The behaviour of S^ Z )(Q)« Q 2 sets in 
only for very small wavevectors: for the largest part of the Brillouin zone excepting 
the neighborhood of Q = n it is rather close to Q 3 ^ 2 . 



Figure 5 

The static structure factors and S^ zz ^ for all chain lengths as a function of 

momentum. The symbol sizes are the same as in figs.l and 6. Crosses stand for 
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V. COMPARISON WITH NEUTRON SCATTERING EXPERIMENTS 

The two magnetic modes in the neighborhood of Q = n have been observed 
by INS in ref. [17-20]. The values of the gaps lead to J « 3.8meV and D/J=0.18. 
These values are very close to those extracted from magnetization measurements 
[36]. The dispersion has been studied throughout the whole Brillouin zone [22]. A 
good fit is obtained through the following formula: 



The parameters are A+ = 2.5 meV, A_ = 1.2 meV, v=9.6 meV and A=6.1 meV. 
Scaling by J =4 meV we plot w±(Q) as continuous lines in fig. 6. For clarity we have 
also plotted our Lanczos points. The edges of two-particle continuum extracted 
from the fits (5.1) are also plotted in fig. 6. The lower dashed line corresponds to 
the (S z = +1) + (S z = —1) continuum while the upper dashed line corresponds to 



Lanczos results compared to experimental fits. The two continuous curves are taken 
from Eq.(5.1). They reproduce the findings of ref. 22. The two dash-dotted lines are 
the two-magnon continuum boundaries obtained also from (5.1). The size of the 
symbols is chosen as in fig.l. 




(5.1) 



(S z = +1) + (S z = 0). 



Figure 6 
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There is very good agreement with the theoretical results in the neighborhood 
of Q=7r where the IP and OP modes are separated by the experimental resolution 
[17-20,22]. In the region n/2 < Q < n the fits (5.1) are also in agreement with 
the theoretical results. The two modes are no longer separated since they are very 
close to each other. At the top of the dispersion curve Q=7r/2 we estimate from the 
values quoted in sect. IV the gaps to be 2.75 J for the IP mode and 2.65 J for the 
OP. This means that the IP mode is at ~ 11 meV and the OP mode at ~ 10.5 meV 
using the value of J quoted above. The crossing (or near crossing) of the two modes 
for Q ps 0.757T precludes their separation in this region of Q. In the regime Q < tv/2 
the separation of the modes increases but the intensity of the scattering decreases 
strongly as is seen from fig. 5. The continuum will appear for Q < 0.3tt and is not 
seen in present experiments [22]. We see it for Q=0.257r and below but there the 
structure factor will be very small. For the whole interval 0.37T — 1.07T our data 
are well reproduced by long-lived excitations as is seen experimentally [22]. The 
asymmetry with respect to n/2 of the spectrum (figsl,6) demonstrates the absence 
of broken translational symmetry. The trend of the magnetic intensity versus Q in 
fig. 5 is that found by INS [22]. Absence of data for Q very close to it forbids us to 
check the square-root Lorentzian behaviour expected for each of the static functions 
and S^ zz \ 

With respect to what has already been done, it would be very interesting to 
separate the two modes near 7r/2 where intensity is not dramatically weak. The 
observation of the two-particle continuum on the contrary should be quite difficult 
since it appears only for Q < 0.257T where magnetic scattering is very weak. The 
observation of very long-lived modes [22] is thus in very good agreement with the 
physics of an anisotropic spin 1 chain. 

Finally we comment on the uncertainties in the numerical data. For the gap 
values there is very good convergence for Q=7r and there one can use all the chain 
lengths at our disposal. For other values of Q it is only at the top of the spectrum 
that an extrapolation can be performed. Concerning the values of Q at which some- 
thing interesting happens (i.e. crossing of the two modes or entering the continuum) 
without being able to make refined extrapolations we can only observe that they 
do not depend much on the chain length. The corresponding numbers should be 
regarded as tentative (i.e. 0.75/t for the crossing and 0.37T for the continuum). 

High resolution INS experiments [30] have revealed that the IP mode is split by 
in-plane anisotropy E J2 i [(Sf) 2 — (Sf) 2 ]. Our results neglecting this further splitting 
thus will apply to the description of actual experiments as long as the resolution is 
not very high. This effect is not expected to change much the numerical figures. 
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VI. CONCLUSION 



We have studied the dynamical properties of a spin 1 chain with single-ion easy- 
plane anisotropy thought to be relevant to the magnetic behaviour of the compound 
NENP. The dynamical structure factor g* 00 ) (Q 5 to) has been computed by a Lanczos 
method on chains of lengths up to 16 spins while some parts of the spectrum have 
been obtained on a 18 spins chain. We have performed a variational calculation to 
obtain the spectrum of elementary excitations in the S z = sector. Both ab initio 
and approximate methods agree very well, elucidating the role of the spin-zero 
defects in the Haldane phase. 

We find that the magnetic excitations are described by two distinct long-lived 
modes throughout most of the Brillouin zone 0.3 — 1.07T where the lowest excited 
state in each sector (S z = and S z = ±1) carries almost all of the spectral weight: 
a single-mode approximation is thus adequate. We find that these modes merge into 
a two-particle continuum for Q < 0.3tt as seen in the frequency dependence of the 
structure factor. We have shown that these results reproduce the INS experiments 
on NENP. The dispersion of the magnetic modes as well as their intensity is close 
to that found theoretically. 

Our results show that in principle it should be possible to observe the crossing 
of the two IP and OP magnetic modes. At the top of the dispersion curve their 
separation and intensity should allow observation. On the contrary the edge of 
the two-particle continuum appear only for very small Q< 0.37T where magnetic 
intensity is very weak as found in recent experiments [22] . We mention finally that 
it would be very interesting to measure the static form factors of an anisotropic 
chain by a method that allows to reach the neighborhood of Q=7r since here we 
are limited by the coarse-graining due to our small chains. Quantum Monte-Carlo 
calculation that works efficiently for static quantities can reach this goal. 
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